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Abstract 



Using a data sample of 21.3 fb~^ collected by the Belle detector, three 

body charmless decays B'^ — > K'^h^h~ have been studied. With no as- 

Jlf^ \ sumptions on the intermediate mechanisms, the following branching frac- 

r~^ ' tions have been measured for the first time B{B^ — > K^t^^t:^) = 

O ■ (58.5±7.1±8.8)xl0-^ and i3(5+ ^ K+K-K+) = (37.0±3.9±4.4) x 10"^ 

We also present the first observations of the decay mode B~^ — > K*^(892)tt^ 

with a branching fraction of i3(B+ -^ i^*0(892)7r+) = (16.7t|^t|Jt||j) x IQ-^ 

O I and the decay mode B^ -^ /o(980)i^"*" with a product branching fraction of 

O,' B{B+ ^ /o(980)E:+) X S(/o(980) ^ vr+vr-) = {ll.Jtl-'jtHtfD x W'^- 

> 
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I. INTRODUCTION 

During the last few years, a considerable amount of new information on charmless 
hadronic decays of B-mesons has been reported, primarily by the CLEO Collaboration. 
The discoveries of the B — *• Ktt and B ^ tttt decay modes H have provided a real basis for 
searches for direct-CP violating effects in the B-meson system. 

However, because of large combinatoric backgrounds, studies of charmless B decays have 
concentrated mainly on two-body decay processes. In this paper, we report the results of a 
study of decays B^ -^ K^h^h^ {h stands for a charged pion or kaon) where no assumptions 
are made about intermediate hadronic resonances. The inclusion of charge conjugate states 
is implicit throughout this report unless explicitly stated otherwise. 

The data sample used for this analysis consists of 21.3 fb~^ taken at the T(4S') (on- 
resonance) and 2.3 fb~^ taken 60 MeV below for continuum studies (off- resonance) . The 
data were collected with the Belle detector |^ operating at the KEKB asymmetric energy 
e"'"e~ collider IHl. 



II. THE BELLE DETECTOR 

Belle is large-solid-angle spectrometer based on a 1.5 Tesla superconducting solenoid 
magnet. Charged particle tracking is provided by a silicon vertex detector (SVD) and a 
cylindrical drift chamber (CDC) that surround the interaction region. The SVD consists 
of three approximately cylindrical layers of double-sided silicon strip detectors; one side of 
each detector measures the z coordinate and the other the r-cj) coordinate. The CDC has 
50 cylindrical layers of anode wires; the inner three layers have instrumented cathodes for 
z coordinate measurements |^. Twenty of the wire layers are inclined at small angles to 
provide small-angle stereo measurements of z coordinates along the particle trajectories. 
The charged particle acceptance covers the laboratory polar angle between 9 = 11° and 150° 
corresponding to about 92% of the full solid angle in the CMS. 

Tracks are fit using an incremental Kalman filtering technique, where individual mea- 
surements found by the CDC pattern recognition algorithm are added successively to update 
the track's parameters and covariance matrix at each measurement surface. This approach 
to track fitting minimizes the effects of multiple Coulomb scattering on the determination of 
the track parameters. Hits from the SVD are associated and included during the last steps of 
this recursion. The momentum resolution is determined from cosmic rays and e^e~ — »• n~^fi~ 
events to be (Jpjpt = (0.30 © 0.19pi)%, where pt is the transverse momentum in GeV/c. 

Charged hadron identification is provided by dE/dx measurements in the CDC, a mo- 
saic of 1188 aerogel Cerenkov counters (ACC), and a barrel-like array of 128 time-of- flight 
scintillation counters (TOF). The dE/dx measurements have a resolution for hadron tracks 
of 6.9% and are useful for ti /K separation for piah < 0.8 GeV/c and piab > 2.5 GeV/c where 
Piab is the laboratory momentum. The TOF system has a time resolution for hadrons that 



is cr ~ 100 ps and provides t[ / K separation for piab < 1.5 GeV/c [Q. The indices of refrac- 
tion of the ACC elements vary with polar angle to match the kinematics of the asymmetric 
energy environment of Belle and cover the range 1.5 < piab < 3.5 GeV/c [0. 

Hadron identification is accomplished by combining the information from these three 
subsystems into a single number using the likelihood method: 

where h stands for the hadron type (vr, K, p). High momentum tagged kaons and pions 
from kinematically selected D*^ -^ D^ti^; D^ -^ K^tt'^ decays are used to determine a 
charged particle identification efficiency of about 90% and a misidentification probability of 
about 8%. 

Electromagnetic showering particles are detected in an array of 8736 CsI(Tl) crystals 
located in the magnetic volume and covering the same solid angle as the charged particle 
tracking system ^. The energy resolution for electromagnetic showers is (Je/E = (1.3 © 
0.07/E© 0.8/^1/^)%, {E in GeV). Neutral pions are detected via their 7r° -^ 77 decay. The 
7r° mass resolution varies slowly with energy, averaging a^n n = 4.9 MeV. 

Electron identification in Belle is based on a combination of dE/dx measurements in the 
CDC, the response of the ACC, and the position, shape and total energy (i.e. E/p) of the 
shower registered in the Csl calorimeter. The electron identification efficiency, determined 
by embedding Monte Carlo tracks in multihadron data, is greater than 92% for tracks with 
Piab > 1.0 GeV/c and the hadron misidentification probability, determined from Ks — > tt+tt" 
decays, is below 0.3%. 

The 1.5 T magnetic field is returned via an iron yoke that is instrumented to detect 
muons and Kl mesons. This detection system, called the KLM, consists of alternating 
layers of charged particle detectors and 4.7 cm thick steel plates. The total steel thickness 
of 65.8 cm plus the material of the inner detector corresponds to 4.7 nuclear interaction 
lengths at normal incidence. The system covers polar angles between 6 = 20° and 155° and 
the overall muon identification efficiency, determined by a track embedding study similar to 
that used in the electron case, is greater than 87% for tracks reconstructed in the CDC with 
Piab > 1.0 GeV/c. The corresponding pion misidentification probability determined from 
Ks -^ TT+TT^ decays is less than 2%. 



III. EVENT SELECTION 

Charged tracks are required to satisfy a set of track quality cuts based on the average 
hit residual and the impact parameters in both the r-0 and r-z planes. We require that the 
transverse track momenta be greater than 100 MeV/c to reduce low momentum combinatoric 
background. All the cuts used for the selection of charged tracks are listed in Table |. 

Charged particles are identified as iiT's or tt's by cutting on the likelihood ratio (PID): 

At large momenta (>2.5 GeV/c) only the ACC and dE/dx are used since here the TOP 
provides no significant separation of kaons and pions. Por all three-body final states except 



TABLE I. Parameters used for selection of charged tracks. 



Parameter Cut value 

Transverse momentum, p^ pt > 0.1 GeV/c 

Impact parameter, AR \^R\ < 0.25 cm 

Impact parameter, AZ \^Z\ < 2.50 cm 

Number of CDC axial hits, Nah Nah > 10 

Number of CDC stereo hits, Nsh Ngh > 5 



the KKK final state the likelihood ratio for kaon candidates is required to be greater than 
0.8. Otherwise the charged track is identified as a pion. For the KKK final states, we 
require PID{K) > 0.5 to maintain high efficiency. 

All charged tracks are also required to satisfy an electron veto requirement, which de- 
mands that the electron likelihood is less than 0.95. In addition, all charged kaon candidates 
are required to satisfy a proton veto: 

The candidate events are identified by using the beam-constrained mass Mbc = 



s/A — P^^ and the energy difference AE = E'^ — A/i/2, where E^ and Pg are the measured 
energy and 3-momentum of the B candidate in the T[4S) rest frame and ^/s is the total 
energy in T(45') rest frame. 

The Mbc signal distribution is well modeled by a single Gaussian function. The central 
value and width of the Gaussian function are determined from the B~^ — »• D^tt"*", D^ -^ K^7i~ 
signal to be 5.2805 GeV/c^ and 2.75 MeV/c^ respectively. The Mbc resolution is primarily 
due to the energy spread of the e"*" and e~ beams and is found to be independent of the 
particular three charged track final state. We use the ARGUS function p to describe 
the background in the Mbc distribution: /{Mb) = vT^-a?exp[— ^(1 — x^)], where x = 
Mb/ El^^^ and ^ is a free fit parameter. 

The AE signal shape is parameterized by the sum of two Gaussians with the same mean. 
The shape of the background in the AE distribution from non-resonant e^e~ —* QQ {q = 
u, d, s, c) continuum events is parameterized by a first order polynomial. The AE shape due 
to background from other B and B decay processes has a substantial dependence on the 
final state being studied. Unless explicitly stated otherwise, the contributions from these 
backgrounds are also parameterized by a first order polynomial. 

In the following, we refer to the "S signal region;" this is defined as: 

5.272 < Mbc < 5.289 GeV/c^; |A^| < 40 MeV. 

IV. BACKGROUND SUPPRESSION 

An important issue for this analysis is the suppression of the large combinatoric back- 
ground that is dominated by qq continuum events. To suppress this background, we use a 
set of variables that characterize the event topology. 
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Since the two B mesons produced from the T(4S') decay are nearly at rest in the CMS 
frame, the angles of the decay products of two -B's are uncorrelated and the events tend 
to be spherical. In contrast, hadrons from continuum qq events tend to exhibit a two-jet 
structure. Figure ||a shows distributions of | cos(^Thr)|, where Oxhr is the angle between the 
thrust axis of the B candidate and that of the rest of the event. The distribution is strongly 
peaked near | cos{dThr)\ — 1.0 for qq events while it is nearly flat for BB events. We require 
I cos{9Thr)\ < 0.80 for all modes under consideration; this eliminates 83% of the continuum 
background and retains 79% of the signal events. 

After the imposition of the cos{6Thr), QQ and BB requirements, the remaining events still 
have some differences in topology that are exploited for further continuum suppression. We 
construct a "Virtual Calorimeter" f^ by dividing the space around the candidate thrust axis 
into nine polar angle intervals of 10° each; the i-th interval covers angles from {i - l)xlO° 
to 2 X 10°. We define the momentum flows, Xi{i = 1,9), into the i-th interval as a scalar 
sum of the momenta of all charged tracks and neutral showers directed in that interval. The 
momentum flow in corresponding forward and backward intervals are combined. 

Angular momentum conservation provides some additional discrimination between BB 
and continuum qq events. In qq production, the direction of the candidate thrust axis, 
Ot, with respect to the beam axis in the cms frame tends to maintain the 1 + cos^(^r) 
distribution of the primary quarks. The direction of the B candidate thrust axis for BB 
events is uniform. The B candidate direction, 6b, with respect to the beam axis exhibits a 
siry^ {9b) distribution for BB events and is uniform for qq events. 



A Fisher discriminant lOl is formed from 11 variables: the nine variables of the "Virtual 



Calorimeter", | cos{9t)\, and | cos{9b)\- The discriminant JF is the linear combination 
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of the input variables, Xj, that maximizes the separation between signal and background. 
The coefficients «« are determined from the Monte Carlo simulation using a large set of 
continuum events and signal events modeled as B~^ —^ K^7i~^7t~. We use the same set of 
coefficients a, for all modes under study. Figure |I]b shows the JF distributions for the Monte 
Carlo signal in the mode B~^ — > K~^7t~^tt~, and the data signal in the mode B~^ — ^ D^ir^ 
followed by D^ —>■ K~^tt^. The JF distributions for Monte Carlo background and below- 
threshold background data for modes comprising three charged tracks are also presented in 
Fig. |I|b. The JF distributions for both the signal and background are fitted to Gaussian func- 
tions. The separation between the mean values of the signal and background distributions 
is approximately 1.3 times the signal width. 

For the Kirn and KKtt final states we make the requirement on the Fisher discriminant 
variable JF > 0.8; this rejects 90% of continuum background with about 54% efficiency for 
the signal. In case of the KKK final states, the continuum background is much smaller and 
we make the looser requirement JF > 0. This rejects 53% of continuum background with 
about 89% efficiency for the signal. 

To determine the dominant sources of background from other 5-meson decay modes we 
use a large set of Monte Carlo generated BB events where both B mesons decay generically. 
Most of the BB related background is found to originate from B^ -^ D^-k^ and B^ -^ 
J /ip{ip{2S))K^ decays. To suppress this type of background we apply the requirements on 



invariant masses of two-particle combinations that are described below. The background 
from the B semileptonic decays is additionally suppressed by the electron veto requirement. 
The most significant background to K'^tt^tt' final state from the B rare decays is found to 
originate from the B~^ -^ rj'K^ followed by rj' —^ p°7. We expect about 3% of the events 
of this type to satisfy all the selection criteria. We find no significant background to the 
K^K^K~ final state from other rare decays of B mesons. 

V. RESULTS OF THE ANALYSIS 

A. B+ -^ K+TT+TT- 

For the K^tt^ti^ final state, we select B candidates formed from three charged tracks 
where one track is positively identified as a kaon and the other two tracks are consistent with 
the pion hypothesis. The resulting two dimensional AE versus Mbc plot for all selected 
K^7^^7r~ combinations is presented in Fig. ^ where the B signal region is inside the box. 
Figure ^ shows the Dalitz plot for candidates in the B signal region. Large contributions 
from B+ -^ D^n+ where D° -^ K+n~- and B+ -^ J/i/j{i/j{2S))K+ where J/^(V'(25)) -^ 
fi~^fi~ are apparent in the Dalitz plot. Modes with J /ip{;ip{2S)) contribute to this final state 
due to the muon-pion misidentification. The contribution from the J/ip{tp{2S)) —>■ e+e~ 
submode is found to be negligible (less than 0.5%) after the electron veto requirement. For 
further analysis we exclude l)° and J /ip{ip{2S)) signals by imposing requirements on the 
invariant masses of two intermediate particles: 

\M{K+'K-) - 1.8651 > 0.100 GeV/c2; 

\M{h+h-) - 3.0971 > 0.070 GeV/c^; \M{h+h-) - 3.686| > 0.050 GeV/c^ 

where h^ and h~ are pion candidates. For the J/ilj{ip{2S)) rejection, we use the muon mass 
hypothesis for charged tracks to calculate M{h^h~). 

The AE and Mbc distributions for the events remaining after the exclusion of these large 
signals are presented in Fig. ^ and Fig. ^d, respectively. Here a significant enhancement in 
the B signal region is still observed; the results of a fit to this AE distribution are presented 



in Table |^. The expected AE and Mbc distributions, which are the sum of luminosity- 
scaled off-resonance data and BB Monte Carlo, are shown as open histograms in Figs. §a 
and ^D, respectively; the contributions from only the BB Monte Carlo sample are shown as 
hatched histograms. In the AE spectrum, the shape of the BB background component is 
approximated as an exponential function with a parameter determined from the BB Monte 
Carlo. As can be seen from the hatched histograms in Figs. ^ and ^d, there is no significant 
contribution to the signal from BB generic decays after the large known backgrounds have 
been removed. 

To examine possible intermediate two-body states in the observed B^ -^ K^-k^-k^ signal, 
we analyze the K^i[~ and vr+vr" invariant mass spectra shown in Figs, ^a and ^, respec- 
tively. To suppress the feed-across between the tt+tt^ and i^+vr" states we require the K^n' 



+ -Tr-^ 



TT ' TT 



invariant mass to be larger than 2.0(1.5) GeV/c^ when making the tt+tt^ [K^ti^) 
projection. The hatched histograms shown in Figs. ^ and ^a are the h^h~ invariant mass 
spectra for background events in the AE sidebands: 
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5.272 < Mbc < 5.289 GeV/c^ and 
-0.080 < A^ < -0.050 or 0.050 < A^ < 0.150 GeV 



scaled by area. 

The K^TT~ invariant mass spectrum is characterized by a narrow peak around 0.9 GeV/c^ 
which is identified as i^*''(892) and a broad enhancement above 1.0 GeV/c^ which is subse- 
quently referred to as i^x(1400). 

In the TT+vr" invariant mass spectrum two distinct structures in the low mass region 
are observed. One is slightly below 1.0 GeV/c^ which is identified as /o(980) and the 
other between 1.0 GeV/c^ and 1.5 GeV/c^ and referred to as /x(1300). The peak around 
3.4 GeV/c^ is consistent with the process B^ — >• XcqK^ i Xco ~^ tt+tt", and is the subject of 
a different analysis |jll|. In this paper we exclude the Xco candidates from the analysis of 
two-body final states by applying the requirement on the tt+tt" invariant mass: |M(7r+7r^) — 
3.4151 > 0.050 GeV/c2. 

For further analysis we subdivide the Dalitz plot area into seven non-overlapping regions 
as defined in Table ||. Regions from I to V are arranged to contain the major part of 
the signal from the B+ -^ ir*°(892)7r+, B+ -^ A"x(1400)7r+, B+ -^ p^{770)K+, B+ -^ 
/o(980)ii''^, and B^ -^ /x(1300)if^ final states, respectively. The area in the Dalitz plot 
where Ktt and vrvr resonances overlap is covered by the region VI, and region VII covers 
the rest of the Dalitz plot. The /S.E and Mbc distributions for each region are shown in 
Fig. 1^ and the results of the fits are summarized in Table |I|. As can be seen from Fig. ^ and 
Table |IT|, the contribution from region VII to the total number of signal events is negligibly 
small. 

TABLE II. Results of the fit to the IS.E distribution for different regions in the i^+vr+vr" 
Dalitz plot. Columns list the definition of each region, reconstruction efficiency from Monte Carlo 
simulation, signal yield and statistical significance. 

Region Mass range, GeV/c^ Efficiency, % Yield, event Significance, a 

I 0.82 < M(E:7r) < 1.00; Af(7r7r) > 1.50 20.7 ± 3.7 SS.ltgJa 6.1 

II 1.00 < M{Ktt) < 1.75; M{inT) > 1.50 19.2 ± 1.3 37.8l|^^ 4.7 

III 0.62 < M(7r7r) < 0.90; M(K7r) > 2.00 16.7 ± 2.4 6.04l|^^ 1.3 

IV 0.90 < MiiTTT) < 1.06; M{Ktt) > 2.00 19.9 ± 3.2 32.0l^:||^ 6.9 

V 1.06 < MiiTTT) < 1.50; M{Ktt) > 2.00 19.6 ± 1.7 25.4t^j^ 4.3 

VI M(E:7r) < 2.00; M(7r7r) < 1.50 14.7 ± 3.3 12.0l|^^ 2.9 

VII M(ii:7r) > 2.00; M(7r7r) > 1.50 16.1 ±0.5 6.451^;^^ 0.9 



B. B+ -^ K+K+K- 

For the selection of -B ^ K^K'^K~ events, we use events with three charged tracks that 
are positively identified as kaons. To suppress the background caused by it/ K misidentifi- 
cation, we exclude candidates if the invariant mass of any pair of oppositely charged tracks 
from the B candidate is consistent with the D^ —>■ Kit hypothesis within 12 MeV (~ 20"), 
independently of the PID information. 
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Figure ^a shows the two-dimensional AE versus Mbc plot for all selected K^K^K~ 
combinations and Fig. ^p shows the Dalitz plot for candidate events in the B signal region. 
Since in this case there are two same-charge kaons, we distinguish the K^K~ combinations 
with smaller, M{K~^K~)min, and larger, M{K~^K~)max, invariant masses. We avoid double 
entries by forming the Dalitz plot as M'^{K^K^)max versus M'^{K^K^)min, as shown in 
Fig. |b. 

The signal from the Cabibbo-suppressed B~^ — > D^pK^, D^p —>■ K^K~ decay mode 
is apparent as a vertical strip in the Fig. ^p Dalitz plot. The notation D^p means that 
the D meson decays to the CP eigenstate. The corresponding Cabibbo-allowed B^ -^ 
DqpTt^, Dqp —>■ K^K~ decays can also contribute to this final state as a result of pion-kaon 
misidentification. The detailed analysis of the decays of type B^ -^ D^pK^ is described in 
ref. [0. 

We exclude candidates consistent with the B^ -^ D^ph^ , D^p -^ K^K hypothesis 
from further analysis by imposing the requirement on the K^K~ invariant mass: 

\M{K+K-) - 1.8651 > 0.025 GeV/c^. 

The AE and Mbc distributions after the exclusion of D mesons are presented in Figs, ^a 
and ^, respectively. A large peak in the B signal region is apparent in both distributions. 
The results of a fit to the AE distribution are presented in Table |TV|. 

The open histograms in Figs. ||a and |b show the M{K^K~)min and M{K^K^)max dis- 
tributions for selected events, respectively; the hatched histograms show the corresponding 
spectra for the AE sidebands: 

5.272 < Mbc < 5.289 GeV/c^ and 
-0.200 <AE< -0.050 or 0.050 < AE < 0.200 GeV, 

scaled by area. The M{K^K^)min spectrum. Fig. §a, is characterized by a narrow peak at 
1.02 GeV/c^ corresponding to the 0(1020) meson and a broad structure around 1.5 GeV/c^, 
which is subsequently referred to as /x(1500). To exclude the possible contribution from the 
B^ -^ XcqK^, Xco -^ K^K^ final state we apply the requirement on the K^K~ invariant 
mass: \M{K^K^) — 3.415| > 0.050 GeV/c^. The study of this final state is described in 
ref. [0. 

For further analysis we subdivide the area of the Dalitz plot into the four non-overlapping 
regions defined in Table |ITl[ Regions I and II are arranged to contain the major part of the 
signal from the B^ -^ (f){1020)K^ and B^ -^ fx{1500)K'^ final states respectively. Regions 
III and IV cover the rest part of the Dalitz plot. The AE and Mbc distributions for each 
region are shown in Fig. ^ and the results of the fit are summarized in Table |ITT|. 

VI. BRANCHING FRACTIONS RESULTS 

To determine branching fractions we normalize our results to the observed B^ -^ D^n, 
D^ -^ K^7i~ signal. Although this introduces a 9.7% systematic error due to the uncertainty 
in the B^ -^ D^tt branching fraction, it removes systematic effects in the particle identifi- 
cation efficiency, charged track reconstruction efficiency and the systematic uncertainty due 
to the cuts on event shape variables. 
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TABLE III. Results of the fit to the A£' distribution for different regions in the K^K^K~ 
Dalitz plot. Columns list the definition of each region, reconstruction efficiency from Monte Carlo 
simulation, signal yield and statistical significance. 

Region Mass range, GeV/c^ Efficiency, % Yield, event Significance, a 

I 1.005 < M{KK)rain < 1-035 24.6 ± 2.5 246lfl ^^ 

II 1.035 < M{KK)min < 2.00 23.3 ± 0.8 84.01^^1 8.9 

III M{KK)min > 2.00; M{KK)max > 3.40 23.9 ± 1.1 15.9l|^^ 3.3 

IV M(KK)™„ > 2.00; M(KK)„ax < 3.40 24.7 ± 0.8 26.1II43 6.7 



We calculate the branching fraction for B meson decay to a particular final state / via 
the relation: 

B{B+ ^ f) = BiB+ -^ D\+) X i3(D° ^ K+n-)^ x ^, (1) 

where Nf and A^dtt are the numbers of reconstructed events for the particular final state / 
and for the reference process, respectively; £/ and Edtt are the corresponding reconstruction 
efficiencies. 

The number of signal events for the normalization processes B^ -^ D^n and D^ -^ 
K~^TT~ is found to be 1137 ± 38 for the K^K^K~ selection requirements and 619 it 29 for 
requirements used for the other Khh combinations with the reconstruction efficiencies of 
31.3% and 16.8% respectively. 



A. Three-body branching fractions 

For branching fraction calculations we use the signal yield extracted from the fit to 
the corresponding Ai? distribution since the Mbc distribution suffers more from BB back- 
ground. The reconstruction efficiency for three-body final states is determined from the 
Monte Carlo simulations of events that are generated to have a uniform distribution over 
the Dalitz plot. 

The branching fraction results are summarized in Table |V|. The first quoted error is 
statistical and the second is systematic. We do not observe a statistically significant signal 
in the K'tt^tt^ , K^K^tt~ or K^K~t:^ final states and place 90% confidence level upper 
limits on their respective branching fractions. 

The following sources of systematic errors are found to be the dominant: 

• uncertainty due to the nonuniformity of the reconstruction efficiency over the Dalitz 
plot. It is found to be 7.6% for K+^-tt+ and 3.7% for K+K+R-; 

• uncertainty in B^ -^ D^tc^ and D^ -^ K^n^ branching fractions: 9.7%; 

• uncertainty in the parameterization of the signal and background shape in AE. This 
is 8.5% for K+Tr-TT+ and 4.7% for K+K+R-; 

• kaon identification efficiency for modes with more than one kaon in the final state: 3% 
per kaon; 
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TABLE IV. Measurement results. Branching fractions and 90% C.L. upper limits for 
B+ -^ K+h+h~ final states. 



Mode 


Efficiency, 


% 


Yield, event 


B,10-^ 


K+7r^7r+ 


16.6 




177 ± 20 


58.5 ±7.1 ±8.8 


K+K+K- 


24.5 




162 ± 16 


37.0 ±3.9 ±4.4 


K^7r+7r+ 


16.2 




3.86l?:i 


<7.7 


K+K+7T- 


13.3 




6.78lt?^ 


<6.0 


K+K-7T+ 


13.2 




oo q+9.49 


<21 



B. Exclusive two-body branching fractions in the i^'+vr+vr final state. 

In the determination of the branching fractions for exclusive two-body final states, we 
take into account the possibility of interference between wide resonances. This requires some 
assumptions about the states that are being observed and, as a consequence, introduces some 
model dependence into the extraction of the exclusive branching fractions. The present level 
of statistics does not permit unambiguous interpretation of the i^x(1400) and /x(1300) 
states and, thus, it is not possible to use the data to fix all of the input model parameters. 
For this analysis we assume that the observed i^x(1400) and /x(1300) states are 0"''+ scalars. 
While this does not contradict the observed signal, some contributions from vector (1^) and 
tensor (2"'") resonances can not be excluded. The uncertainty related to this assumption is 
included in the model-dependent error described below. We ascribe to the Kx{l4:00) state 
the parameters of 7^^(1430) (M = 1412 MeV/c^, T = 294 MeV) and to /x(1300) state the 
parameters of /o(1370) (M = 1370 MeV/c^, T = 400 MeV) [|13l. 

For further analysis we make following assumptions: 

• The observed signal in the i^+vr+vr" final state can be described by some number of 
two-body final states. We restrict ourselves with the following set of exclusive two-body 
final states: ir*0(892)7r+, irx(1400)7r+, p^{770)K+, fo{980)K+ and /x(1300)K+. We 
enumerate these final states as 1 through 5 in the order they are mentioned above. 

• Given this set of two-body final states, we determine the exclusive branching fractions 
neglecting the effects of interference. The uncertainty due to possible interference be- 
tween different intermediate states is included in the final result as a model-dependent 
error. 

In order to extract the signal yield for each two-body final state, we perform a simul- 
taneous likelihood fit to the AE distributions for the seven regions of the K~^it~^7t~ Dalitz 
plot (see Fig. ^ and Table IT]). We express the expected number Uj of signal events in the 
j-th region of the Dalitz plot as a linear combination: 



n 



/ . ^ij^^ii 



'3 /1^^«J^'«5 



where Ni is the total number of signal events in the z-th two-body final state and eij is the 
probability for the z-th final state to contribute to the j-th region in the Dalitz plot. The 
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Eij matrix is determined from a Monte Carlo simulation and includes the reconstruction 
efficiency. This procedure takes correlations between different channels into account when 
determining the statistical errors. 

The results of the fit are summarized in Table 0. Combining all the relevant numbers and 
using Eq. 1, we calculate the product of branching fractions B{B~^ -^ Rh~^) x B{R -^ h^h~), 
where R denotes the two-body intermediate resonant state. The branching fraction result 
for the B~^ — >■ i^*''(892)7r+ final state is in agreement with the results of a separate study of 



B meson decays to the pseudo-scalar and vector final states [Q. 

We present three types of error for the branching fractions: the first error is statistical, the 
second is systematic, and the third reflects the model-dependent uncertainty. In general, 
the model-dependent error is due to uncertainties in the effects of interference between 
different resonant states. We estimate this error by means of a B~^ -^ K^tt^tt~ Monte 
Carlo simulation that includes interference effects between all final states mentioned above. 
We vary the relative phases of resonances and determine the signal yield using the procedure 
described above. The maximal deviations from the central values are used as an estimate 
of the model dependence of the obtained branching fractions. 

TABLE V. Results of the simultaneous fit to the K^tt^tt^ final state. 
Two body mode Efficiency, % Yield, events Significance, a BB+^Rh+ ^ ^R^h+h~-> 10^^ 



K*0(892)^+ 


19.1 


oo c+8.50 
"J°-"J-7.90 


6.2 


-11 -1+2.5+1.4+2.0 
lJ---'^-2. 3-1.4-3. 9 


Kx(1400)^+ 


17.0 


oq 1+10.8 


4.1 


-ir, 7+3.5+1.8+2.9 
-■-^•'-3.4-1.8-5.8 


p0(770)if+ 


18.9 


1 7c:+8.60 
-'-■'^-7.38 


0.2 


<9.6 


/o(980)K+ 


19.3 


m.^tlf. 


6.0 


-1-1 7+2.5+1.5+4.1 
11- '_2. 7-1. 5-1.0 


/x(1300)K+ 


17.3 


33.6tJ2.J 


3.2 


-in 7+3.9+1.4+6.9 
-'^'-'•'-3.5-1.4-2.8 



We find that the model-dependent errors associated with the wide resonances intro- 
duce significant uncertainties into the branching fraction determination. In the case of the 
p'^{llQ)K^ final state, this effect is enhanced by the smallness of the signal itself. The ef- 
fects of interference of a small p°(770) signal with the large /o(980) signal or with a broad 
/x(1300) resonance could result in as much as a doubling of the observed p°(770)i^"'" final 
state or its total suppression. Since we do not observe a significant signal in this channel, 
we report a 90% confidence level upper limit. 

C. Exclusive two-body branching fractions in K^K^K~ final state. 

In the case of the three charged kaons final state we clearly observe the presence of the 
0(1020) meson plus a very broad /x(1500) structure that we currently cannot interpret 
unambiguously. It could be a complex superposition of several intermediate states and 
some contribution from non-resonant B^ — > K^K^K~ decay is also possible. For our 
study of systematic and model-dependent uncertainties, we construct a simplified model 
and parameterize the /x(1500) structure as a hypothetical scalar state with M = 1500 
MeV/c^ and F = 700 MeV. We find qualitative agreement between the experimental Dalitz 
plot distribution of the signal events and that obtained from the Monte Carlo simulation 
with this simple model. 
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Then we extract the signal yield for the two-body final states: B~^ -^ 0(1020)-^'"'" and 
the so-called B~^ -^ /x(1500)i^"'", which is, in fact, all of the remaining signal. We follow 
the same procedure as we used for the i^+vr+vr" final state. The signal yields are determined 
from a simultaneous fit to the AE distributions for four separate regions of the K^K^K~ 
Dalitz plot (see Fig. |^ and Table [TTID . The results of the fit are summarized in Table |V^. The 
branching fraction result for the B^ -^ 0(1020)7^"'" final state is in good agreement with the 
results of a dedicated analysis of the B+ -^ 0(lO2O)i\:+ and B+ -^ 0(lO2O)/s:*+(892) final 
states PI, B{B+ -^ (/.(1020)i^+) = (10.61?;^ ± 2.2) x 10"^. This latter number should be 
considered as the current "official" Belle result for the B^ —>■ 0(1020)7^"'" branching fraction. 



TABLE VI. Results of the simultaneous fit for the K^ 


-K+K 


final state. 




Two body mode Efficiency, % Yield, events 


Significance, a 


Bb+- 


-+fiA'+ X Br-^k+k- 


, 10-'^ 


0(1O2O)K+ 25.1 30.lt^:i 
/x(1500)K+ 22.1 122t\li 


6.4 
12.1 




a 7n+1.6+0. 8+0.4 
"■"J-l. 5-0. 8-0.4 
on 0+3.8+3.9+1.5 
"-"-'•°-3.7-3.9-1.5 





We determine the model-dependent error the same way as we did for the i^+Tr+vr" final 
state. In the case of the K^K^K^ final state the model dependent error is found to be 
much smaller then in K^tt^tt~ final state. This is mainly due to the narrow width of the 
0(1020) meson. 

VII. DISCUSSION & CONCLUSION 

The high quality Belle's tx/K separation allows us to measure for the first time the 
branching ratios for the three-body modes B{B^ -^ K^tt^tt^) = (58.5 ± 7.1 ± 8.8) x 10~^ 
and B{B+ -^ K+K-K+) = (37.0 ± 3.9 ± 4.4) x 10"*^ without assumptions about particular 
intermediate mechanisms. CLEO |]16| and BaBar fl^ have previously placed upper limits 



on non-resonant three-body decays; the reported numbers for B^ -^ K^ti^ti (CLEO: 
B{B+ -^ K+TT+TT-) < 28 X 10-^ BaBar: B{B+ -^ K+n+n-) < 66 x 10"*^) and B+ -^ 
K+K+K- (CLEO: B{B+ -^ K+K+R-) < 38 x 10"^) are lower than those presented in 
this paper. A comparison of the applied selection criteria shows that CLEO and BaBar 
restricted their analyses to the region of the invariant masses above 2 GeV/c^ for any pair 
of the particles. This requirement effectively removes most of the low mass resonances 
that provide the dominant contribution to our observed signal. They assume a uniform 
distribution of events over the Dalitz plot to obtain the limits quoted above. 

The upper limits reported here for the K~7r^7^^, K^K~tt^ and K^K^t^^ modes are 
considerably more restrictive than previous limits from CLEO [|16| {B{B^ — > K^-k^tt^) < 
56 X 10-^ B{B+ -^ K+K-TI+) < 75 x 10"^) and OPAL ^ {B{B+ -^ K+K+n-) < 
87.9 X 10-^). 

Significant signals are observed for the first time in the decay modes B^ -^ 7^**^(892)71 "*" 
and 5+ -^ fo{980)K+ with branching fraction products of B{B+ -^ K*^{892)n+) x 
BiK*\892) -^ K+n~) = (ILltaJtl^tl^) x lO'^ and S(5+ -^ /o(980)ir+) x S(/o(980) ^ 
TT+TT^) = (11.7J]2!7ti;5li;o) X 10^^, respectively. The latter final state is of interest since this 
is the first observation of a i? decay to a charmless scalar-pseudoscalar final state. The signif- 
icant enhancement in the K~^tt~ invariant mass spectrum above the i^*(892) mass (referred 
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to as ii'x(1400)) agrees with the scalar i^Q(1430) hypothesis. This is also in agreement with 
theoretical predictions 0,^ for the B^ — >■ i^Q(1430)7r+ branching fraction made based on 
the factorization model. Nevertheless, we cannot exclude some contribution from the tensor 
K*(1430) state. 

Large uncertainties arise in the interpretation of the peak with a tt+tt^ invariant mass 
around 1300 MeV/c^ in the Ktttt system. There are two known candidate states: the 
/2(1270) and /o(1370) 0. Attributing the peak to the /o(1370), with its rather small 
coupling to TT+TT" pl|, would lead to an unusually large branching fraction for a charmless 



B decay mode. On the other hand, as recently shown in \72\, factorization model predicts a 
very small branching ratio for B^ -^ f2{1270)K^. If our observation is, in fact, due to the 
/2(1270), this would provide evidence for a significant nonfactorizable contribution. 

We cannot identify the broad structure observed in the B~^ — >■ K^K^K~ final state above 
0(1020) meson. It is hardly compatible with the presence of single scalar state either /o(1370) 
and /o(1500) |TB[. We also cannot exclude the presence of a non-resonant contribution or the 
case of several resonances contributing to the excess in the K^K~ invariant mass spectrum 
seen around 1.5 GeV/c^. 

We find that effects of interference between different two-body intermediate states can 
have a significant influence on the observed two-particle mass spectra and a full amplitude 
analysis of three-body B meson decays is required for a more complete understanding. This 
will be possible with increased statistics. 
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FIG. 9. The AE (left) and Mbc (right) distributions for different regions of the K^K^K 
DaUtz plot. The plots from the top to bottom correspond to regions I to IV, respectively. 
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